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\Q . Abstract 

Standard methods of nonlinear dynamics are used to investigate the stability 
^ | of particles, branes and D-branes of abelian Born-Infeld theory. In particular 

the equation of small fluctuations about the .D-brane is derived and converted 
into a modified Mathieu equation and - complementing earlier low-energy inves- 
tigations in the case of the dilaton-axion system - studied in the high-energy 
domain. Explicit expressions are derived for the S-matrix and absorption and 
reflection amplitudes of the scalar fluctuation in the presence of the D-brane. The 



O 



results confirm physical expectations and numerical studies of others. With the 
derivation and use of the (hitherto practically unknown) high energy expansion 
of the Floquet exponent our considerations also close a gap in earlier treatments 
D ! of the Mathieu equation. 

> 

x. 

c3 ; 1 Introduction 

Recently Born-Infeld gauge theory has attracted considerable interest as the 
bosonic light-brane approximation or limit of superstring theory]!], and has 
turned out to be a simple and transparent model in this context 0. Branes, 
defined as extended objects in spacetime, can be fundamental or solitonic. The 
connection of these branes with a U(l) gauge field was motivated by the pres- 
ence of this field in the massless part of the spectrum of open strings, and by 
realising that branes with open strings attached to them which satisfy Dirichlet 
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boundary conditions, or more generally one brane attached to another, can be- 
come classically stable, solitonic objects. It is for this reason that the dynamics 
of D-branes in Born-Infeld theory is being studied in detail and generalised 
01) 0)0)0; HI- Since, in general, a brane may or may not be a solitonic configu- 
ration or BPS state, the exploration of this question deserves particular attention. 
It is often stated that a brane is BPS in view of the vanishing of a fraction of 
the supersymmetry variation of the associated gaugino field. However, since BPS 
states (as classically and topologically stable states) and Bogomol'nyi bounds 
have been studied in great detail in a host of other theories, and the approach in 
these is practically standard, one would like to understand aspects of Born-Infeld 
particles in a similar way, also because it is not absolutely clear that .D-branes 
are solitons of string theory in precisely the same way as more familiar topolog- 
ical solitons in field theory. Therefore our first intention in the following is to 
study Born-Infeld particles with standard methods of nonlinear dynamics in the 
simplest case of a flat spacetime. We begin with the free Born-Infeld particles, 
i.e. Blon and catenoid [|]]. Using a scale transformation argument || we show 
that these static configurations - which differ from ordinary solitons of nonlinear 
theories in requiring a special consideration of source terms or boundary condi- 
tions (cf. also (10|, [[□]]) - require the number of space dimensions p to be larger 
than 2. We assume spherical symmetry and study the local stability of these 
configurations by considering the second variational derivatives of their respec- 
tive actions. Our conditions for stability are a) that the eigenfunctions of the 
corresponding operator be square integrable, and b) that the charge e be fixed, 
with angular fluctuations ignored. We then consider the case of the scalar field 
corresponding to a single transverse coordinate coupled to the gauge field (here 
only the electric component), i.e. the catenoid or brane with associated open 
fundamental string. We distinguish between two types of arguments in deriving 
the linearised fluctuation equation, and infer the stability of this stringy D-brane. 
In ref. |L2| an explicit and detailed consideration of the Bogomol'nyi bound in a 
special model of Born-Infeld theory has been given where the central charge of 
the supersymmetry algebra plays the role of the topological or winding number 
of ordinary solitons. 

Our second intention in the following is the explicit study of the small fluc- 
tuation equation about the D3-brane in the high energy domain. This equation 
with singular potential has the remarkable property of being convertible into a 
modified Mathieu equation which depends only on one coupling parameter which 
is a product of energy and electric charge. The S*-matrix for scattering of the 
fluctuation off the brane can be obtained in explicit form. The Z)3-brane is 
therefore one of the very rare examples allowing a detailed study of its properties 
with explicit expressions for all relevant physical quantities in both low and high 
energy domains. We therefore expect that also 5-duality can be uncovered and 
studied in this case (although we do not attempt this here). Various other Dp- 
brane models have been discovered recently whose small fluctuation equations 



can be reduced to modified Mathieu equations [13], [14], [H| which have then been 



investigated mainly by computational methods. For the AdS/CFT correspon- 
dence the logarithmic corrections to the low energy absorption probability are 
of particular interest, since these permit a direct relation to the discontinuity of 
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the cut in the correlation function of the dual two-dimensional quantum field 
theory. The first such logarithmic correction to the absorption probability was 
originally obtained in refs. |T6|, |TT| , IS| without resorting to the use of Mathieu 
functions. Subsequently the authors of ref. [13| considered the modified Mathieu 
equation and used computational methods to generate explicit series expansions 
up to several orders for the low energy absorption probability. In 
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a different 

choice of expansions was considered to obtain leading expressions more easily. It 
is natural to supplement such investigations by exploring also the high energy 
case, the first such consideration being that of ref.||15|. The analytical high en- 
ergy results obtained in the following and the complementary low energy results 
of ref. 
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(we also demonstrate how the S'-matrices are related) are therefore 
directly applicable to these. Singular potentials have been studied from time to 
time, and have mostly been discarded as pathological. It seems, however, that 
their real significance lies in the context of curved spaces with black-hole type of 



absorption |20[ . 

Sections 2 and 3 deal with the Blon and the catenoid, sections 4 and 5 with 
the Bogomol'nyi limit of the D3-brane and the derivation of the linearised fluc- 
tuation equation about it. In section 6 we consider this equation in detail in the 
high energy domain and calculate the rate of absorption of partial waves of the 
fluctuation field by the brane. That this absorption occurs is attributed to the 
singularity of the potential. The absorptivity part of the paper may be looked at 
as the high energy complement to the low energy case of ref. |19| with the same 
expression of the S"-matrix. All these calculations require a matching of wave 
functions. In the low energy S-wave case simple considerations of Bessel and 
Hankel functions suffice as was shown in refs. [21|. The low energy limit is, in 
fact, independent of the choice of matching point, as was shown recently 
Our considerations here, however, are general. 



2 The Blon 

We consider first purely static cases and write the Lagrangian of the static Blon 
in p + 1 spacetime dimensions (cf. ||) 

£ 

L = f d^xC, C = l- y/l- (<9 4 0) 2 -£ p e0<5(r), S p = ^ (1) 

(i — 1, • • -,p) with the charge e held fixed by the constraint 

1 f , did) , , 

e + — dai i g = 2 
Z P J - (9,0)2 

Eq. (D is the La grangian one obtains from the world brane action of the pure 
Born-Infeld U(l) electromagnetic action reduced to the purely electric case with 
field Ei = doAi — diAo and no transverse coordinate. The field is assumed to 
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depend on the world brane coordinates x M , fi = 0, • • -,p. The static Blon equation 
of motion is 

di\ , ^ , i =-V->'(r) (3) 



V 1 - w) s 



In the special case p = 3 the classical SO(3) symmetric solution, called a Blon, 
is given by 



V e 



«/»c(0) - r + 



~,5 



10e 2 



(4) 



and C (O) = |-B(|,|).e5 = 1.854074677.ei, B being the Bernoulli function. It 
is easily verified that this solution satisfies the constraint (0) for any value of r. 
Defining E = — V0 C (so that <ft c = A with dA (xi,t)/dt = in the static case), 
and defining D = ||j = -/=p we have (with F j = 

<9£ 1 
T 00 = F 0i — - - C = E ■ D - C = - 1 + 4vre05(r) (5) 

oF 0i VI - E 2 

The energy H c of the Blon (obtained by integration over R 3 ) is then found to be 
finite, i.e. 

H c = J dxT 00 = 47r(3.09112).e5 (6) 



V 

-1 



and in p dimensions the total energy of the Blon scales correspondingly as 
The finiteness of the energy depends on the minus sign in (|l|) and so with (|3]) on 
the relation 



,2 J2 



i - m 2 = = —^=t (7) 



e *Jl + r 



for < r < oo. It may be noted that by defining D such that the left hand side 
of eq.@ is diD^ the singularity of the right hand side is associated with D rather 
than with E which is the decisive difference between Maxwell and Born-Infeld 
electrodynamics. A similar observation applies to the catenoid equation below. 
The energy of the Blon is seen to be independent of its position which hints at 
the existence of some kind of collective coordinate. However, exploring this point 
further is expected to be difficult since a moving charge generates a magnetic 
field, and hence the electric field alone would not suffice. 

We can use a scaling argument || to show that here finite energy configu- 
rations require p to be larger than or equal to 3. Under a scale transformation 
x — > x 1 = Ax, 0(x) — > 4>\(x) = 0(Ax), di<j)(x) — > [di4>(x)}\ = A<9j0(Ax). The charge 
e defined by the constraint (0) also changes under the scale transformation, i.e. 

1 r deride 



X p ~ 2 Z p J J i _ A 2 (<9j0) 2 
In particular for p = 3 and radial symmetry 



et 3) = g 1 1 0) 
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and for arbitrary values of A the r-dependence drops out only if the limit r — > oo 
is taken in the evaluation of the integral. Then 

>=3) i 

r ^ i (10) 
e A 

But also e\=i = e for any r. If C is stable and 7^ 0, the energy must be 
stationary for A = 1, i.e. (dH c /d\)x=i = 0. From this one finds that p > 3. 
Also (d 2 H c /d\ 2 )x=i > for p > 3. Eqs.(||) and (|T0D show that changing the 
scale changes both the charge and the energy, i.e. if the charge were variable, 
one could lower the energy and hence the configuration could be unstable. But 
fixing the charge (e.g. by a quantisation condition) no instability is implied by 
the scaling condition. 

We investigate the stability of the Blon further in the special and exemplary 
case of p = 3 by considering the second functional variation of the static La- 
grangian evaluated at c ( r )- This can be written and simplified in the following 
form (ignoring total divergences on the way) 

5 2 L = - f d 3 x5<f)A5<j), (11) 
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where 

A 



- d I d . - d W fl. 



Id r 2 d 
The operator A can also be written 



(12) 



a f 1 d 2 d 6 ,,2 d \ , , 

= 0/2)3/2 \^d^ r fo~r c ~fo) ( ' 

The classical stability of </> c is therefore decided by the spectrum {u n } of the small 
fluctuation equation 

-hi(i-^i^ = ^ n (14) 

We explore first the existence of a zero mode ipo, i.e. the case u — 0. In this case 

r 2 d 



(1 - 0^) 3 /2 dr 

and so with ^0(00) = 



Vo = C (15) 



C [°° x 4 d [°° dx d(j) c 

^ = -7*1 (IT4j^ rfx = ~^X (TT4j^ = "^ (16) 
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The derivative of the classical configuration C with respect to the charge e indi- 
cates that a perturbation along d(j) c /de around c leaves the static action invariant, 
i.e. <p c (e,r) and c (e + 6e,r) have the same action since 



<90 c (e + Se, r) 



d(5e) 



__ d(p c (e,r) 
Se=o de 



We now show that the operator A does not possess negative eigenvalues, and 
that therefore the Blon is a classically stable configuration. We let ip n be an 
eigenfunction of the operator A. Then 



r „ r°° 
J d xipnAtpn = -Art j drip n 



d r dip n 
~dr (1 - 0' c 2 )3/2 dr 



oo 



-4*1 dr{U, ^ ^ ^ 



dr rn (1 - 0' c 2 )3/2 dr (1 - 0' c 2 )3/2 ^ dr 



OO 



F+ n *<r«*Fx£) (17) 



where F := F(r)|g° and 



The second term on the right hand side of eq.(|TT|) is strictly positive. Hence non- 
positive eigenvalues imply a nonvanishing negative value of F. /,From the condi- 
tion / °° ip^r 2 dr < oo, (i.e. ip n ~ l/(r 1+e ),e > 0) it follows that r 2 ip n dtfj n /dr — > 
with r — ► oo, so that 



Fir) ~ — 47rr -0 n — > 

dr 



and F(oo) = 0. Hence 

F = -F(0) ~ 4vre 



3 V>n 



r 4 c?r 

As r — > eq.(ITJ) becomes 



(19) 



_ L^LL^ = ^ (20) 

r 2 dr r 4 dr n e 3 n 

In the case of the zero mode 

~ d + C 2 r 5 , r -> (21) 



In this case F = 207re 3 CiC2. For C1C2 < this is in full compliance with (|T 
and (|j) from which we obtain 



8eV2 5e 3 J ' 
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For oo n 7^ the small-r behaviour of ifj n is 

1 uj. 



C I 1 - ^ + 0(r- 



(22) 



so that 



4 



1 7 







r=0 



Thus the conclusion is that for all eigenfunctions ip n 



(23) 



This inequality excludes the possibility of the existence of negative eigenvalues. 
Hence the Blon is in this sense classically stable. 



3 The catenoid 

The Lagrangian of the static catenoid in p + 1 spacetime dimensions and with a 
source term is given by (cf. H) 

L = J d p xC, C = l- /l + (<%) 2 - Eprg- V(r) (24) 

where the signs have been chosen such that the energy is positive. Here the scalar 
field y(xi,t) originates from gauge field components A a for a = p + 1, • • (d — 
1), d =dimension, which represent transverse displacements of the brane; here we 
consider the case of only one such transverse coordinate, i.e. y, all d — p — 1 of 
which are essentially Kaluza-Klein remnants of the d = 10 dimensional N = 1 
electrodynamics after dimensional reduction to p + 1 dimensions. The Euler- 
Lagrange equation of the static catenoid y c (static meaning dy(xi,t)/dt = 0) is 
given by 



diVc 



di I , " 1 » c = I = Eprg-tyr) (25) 



,2 



so that after integration 



or for r > r 



'1 + (?/c 
1 + 



-rv (26) 



p-i . — - r p— i 



^w ^P ' ^ 1+ ^ =(+) 7^^F (27) 

In the case of the catenoid without source term the right hand side of eq.([26|) 
can be taken to originate from a boundary condition such as V ■ = 0. The 
domain r < ro is the nonsingular throat region (i.e. y c {ro)) is finite). One may 



observe that the singularity on the right hand side of eq.fl25D is associated with 
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the entire expression on the left whereas, like <9j0 c in the Blon case, so now here 
Vy c is finite, i.e. the p-brane or single throat solution is given by 

_p-i 



Vc(r) =(+) 



dr 



,2p-2 



P 2P-S 




(2? 



Thus y is double valued. The two possible signs can be taken to define a brane 
and its antibrane. We show at the end of this section that the solution with 
the minus sign is the minimum of the action and the solution with the plus 
sign the maximum of the action. This function is finite at r = ro and can be 
expressed in terms of elliptic integrals. For ro = 1 it is even simpler and has 

the value y c (X) =(+) 7^^(^) where /C is the complete elliptic integral of the 
first kind. Plotted as a function of r, y c {r) is a monotonically decreasing function 
starting from ro; pictured on a 2-dimensional space it looks like an inverted funnel 
(i.e. the surface swept out by a catenary with boundaries at the openings), thus 
suggesting the name catenoid. As pointed out in ref.0, the two possible signs 
of the square root allow a smooth joining of one such funnel-shaped branch to 
an inverted one connected by a throat of finite thickness, the resulting structure 
then representing a brane-antibrane pair. This brane-antibrane pair is joined by 
the throat of finite thickness ro and finite length. In fact, we can rewrite eq.(|28|) 
in terms of y c (x) = y c (rox),x = — , and for the special case of p = 3 as 



y c {x) 




where x > 1 and we used formulae of ref . [23 
the periodic function 



(29) 



Inverting this expression we obtain 



x{y) 



cn 



(30) 



Plotting this expression with x as ordinate, one obtains the picture of a cross 
section through a chain of periodically recurring funnel-shaped structures to the 
one side of the throat, i.e. the series U U U U ■ ■ •• representing a series of brane- 
antibrane pairs along the abscissa. 

Proceeding as in the above case of the static Blon and calculating the second 
variational derivative we obtain 

1 
2 



6 2 L 



d p x5yB5y 



(31) 



where for r > ro 



B 



dj- 



[i- 

1 d 



(dty, 



\2lV2 



di - di 



diy c djy c 



[1 + (d lVc 



\ 21 3 / 2 



d 



r 2 dr 



,2\ 3/2 dr 

tic I 



+ 



1 d 



r*fl 2 d 



r 2 ^ r 



dr 



(32) 
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The operator B can also be written 

p 1 f j_ d_ 2 _d 6 /2 d_\ 

(1 + y' 2 fl 2 \ r 2 d/ dr + r^ c dr J 1 j 

Since the gauge field components A a , a = p + 1, • • •, d — 1 (of which we retain only 
one), are dynamical, the Lagrangian in the nonstatic case is 



£ = 1 - y/l-idpyWy) - Eprg-^r) (34) 

and we can obtain the same condition of stability by considering the dynamical 
fluctuation rj, i.e. 



y(*,x) =y e (r)+r}(t,x), rj = i(r)e 



and linearising the time-dependent Euler-Lagrange equation. The square in- 
tegrable perturbations £(r) are the socalled "L 2 deformations" of ref.Q. The 
classical stability of y c is therefore decided by the spectrum {u} of the small 
fluctuation equation 

1 d r 2 did f - (r 4 - r 4 ) 3 / 2 1 d 

r 2 dr(l + ^)3/ 2 d/~r 2 dr\ t+j r 4 } d^ ~ ^ {6b) 

We explore first the existence of a zero mode ^o, i.e. the case uj = 0. In this case 

r 2 d 



(1 + y ,2 fl 2 dr 
and so in the case p = 3 and r > r$ 



^o = C (36) 



ipo = C dx— j-^ = — — 37 

so that 

or z 

Here again the derivative of the classical configuration y c with respect to the 
parameter tq 2 is indicative of stationarity of the action in a shift of rg. 

We now demonstrate that the operator B with the minus sign has no negative 
eigenvalues, and that therefore the free catenoid is a classically stable configura- 
tion like the Blon for fixed throat radius vq. Then 



p A /*oo 

/ d 3 xipBip = An l drip 

J Jrn 



d r dip 
~dr (1 + y' 2 f/ 2 ~dr~ 



r 2 dip 00 f°° r 2 / dip\ 

(l + y' c 2 ) 3 / 2 ^ ro ~ 47r X (l + ^ 2 ) 3 / 2 



2 



(4 4)3/2 ^ 



oo + /-oo (V 4 _ ro 2 ") 3 / 2 



r 4 — r ) ' / dip 



2 



(-) An dr- . 

ro Jro T \dr 



(3? 
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where we used eg. (JZTD . The second term is always positive if the upper sign is 
chosen. The first term vanishes at infinity with / drr 2 ip 2 < oo, since 

(r 4 — r 4 ) 3 / 2 dib r ~,oo 9 dib 

-4tt^ -p— — ib—^- ~ -47rrV-r- ^ °- 

r 4 dr ar 

On the other hand, in the case r — > r , we have 

( r 4 4\3/2 ^ ^ 



As r — > r eq.(p^) becomes 



8 rf (r-roP?=^ (39) 



3/ 2 <ir <ir 



In the case of the zero mode ibo with uj = the considerations are analogous to 
those of the Blon case and the sum of the two terms in eq.fl3"8|) vanishes. In the 
case of uj ^ we therefore have 



and 

lim (r - r ) 3/2 ip— = --C 2 . lim r 3/2 (r - r ) = 
r- >ro dr 8 r ^ r o 

This proves that for all eigenf unctions ib 

< if>\B\*f> > > 0. 

Thus B has no negative eigenvalues, and the free throat is classically stable with 
fixed ro for the sign chosen as in eq.(|38D. Obviously the operator B with the 
plus sign has no positive eigenvalues, which means that we have the maximum 
of the action. Of course, if we change tq (and so consider a different theory), 
the expectation value of B also changes. One should note that the free throat 
we discuss here is that with vanishing gauge field. The double valuedness of the 
solution of eq . (p5|) implies that if one solution is classically stable, the other one 
is not. Thus a multi-throat solution constructed from these by matching both 
solutions, if it exists, like the brane-antibrane solution of ref.[[|, is expected to 
be unstable in view of negative as well as positive eigenvalues, and is therefore 
neither a maximum nor a minimum of the action. In fact, as argued in ref.0] 
(after eq.(132)) equilibrium between these should not be possible. The reason for 
this is that a symmetrical configuration, symmetrical about the plane 23 = for 
instance, implies d^y = there. Evaluating the stress tensor element T33 (even 
for vanishing gauge field), one obtains a negative quantity which is interpreted 
as implying an attractive force between the brane and its antibrane in this sym- 
metrically constructed configuration. This is, in fact, the general instability of 
this configuration discussed in ref.pj. 
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4 Coupled fields: The D— brane in the Bogo- 
mol'nyi limit 

In the case of coupled fields and y (the former with source, the latter without), 
the Lagrangian of the static case is (cf.|§]) 



L = J d p xC, C = l-Q- £ p e0<5(r) 



Q 



1 - {d^f + {d iy f + {d^Ayf - W) 2 (%) 5 



(40) 



^From the first variation of L we obtain the coupled equations of the fields and 
y, i.e. from 

SL= - j5<l ) d~[d i <p-(V<P-Vy)d l y + (Vy)%<P}d p x 
Syd~[d iy + (V0 • Vy)d4 - (V<P) 2 d iy }d p x 



- S p e J 5<f)5(r)d p x (41) 

(ignoring total divergences). 

The source term of the electric field again suggests spherical symmetry. In 
deriving the two coupled Euler-Lagrange equations one new constant c (apart 
from e) arises in the integration of the catenoid equation, i.e. 

d r (r^J^) = 0, r^-^- - ~ 



d(d r y) J d(d r y) 

We have no source term of the y field because, as before, the appropriate effect 
is provided by the boundary condition defining the width of the throat. The two 
equations with spherical symmetry are found to be 



<P' e -y' 



so that 



Then 



_ e = 1 
y' c a 



(42) 



(43) 



a 2 



'■■ ■ < I; 2 5 (V ) r 2(p-l) , 1 2 

- — 2 V 1 — a 1 - — 2 V 1 — a 

Thus the family of solutions can be parametrised in terms of the single parameter 
a as already pointed out in ref. 0. This parameter is seen to interpolate between 
the two types of static solutions. The solution y of fl3"5[) for various values of a 2 
is now the p-brane, i.e. 

+ roc \ 

y(r) =(— ) ae / dr = (45) 

Jr v / r 2(p-l) _ r 2(p-l) 
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where f§ p ^ = e 2 (a 2 — 1) and for the solution to make sense we must have a 2 > 1. 
If ae in eq. (|36|) is replaced by — ae, the expression represents the corresponding 
antibrane. Taking e 2 — > 0, a 2 e 2 — > const, the electric field is eliminated and we 
regain the free catenoid solution. In approaching the limit a 2 — > 1 the width of the 
throat becomes infinitesimal with nonvanishing electric field and the configuration 
can then be considered to be a fundamental string, as argued in ref.0. We 
distinguish between three cases: 



\a\ < 1 



\a\ > 1 



dx f°° dx 

--,y = a 



l-a 2 + x 4 /e 2 J r J\-a? + x A /e 



doc f doc 

--,y = ae 



- ■ 

:,2 



A _ 



\Jx A — r 4 ' Jr i/r" 1 — ro H 

±1 : 0=£y = ±f (46) 



We see that for a 2 = 1 eq. (153) becomes the first order Bogomol'nyi equation or 



linearised field equation for y (as in ref. @) 

F 0r ± f r = (47) 

where Fq t = E c is the static electric field. This is the same equation as that 
obtained from the vanishing of the supersymmetry variation of the gaugino field 
E for half the number of 16 supersymmetries (for d — 10 and p = 3) e+, e_ of e 
for which <5£ = 0, i.e. 

= 0, ^ 



where - as discussed in the literature |24| - e is the constant spinor of the su- 
persymmetry variation and e± are its chiral components. Thus a 2 = 1 implies 
BPS configurations, wheras those with a 2 ^ 1 are non-BPS. Taking a 2 = 
in eq. (|36l) we regain the Blon configuration as a local minimum of the energy 
whereas for vanishing electric field one expects a local maximum, i.e. a sphaleron 
configuration (as pointed out in @). 

Next we investigate the second variation of the static L with spherical sym- 
metry. We set 

5 2 L = ^J \5<pM5<p + SyNSy + 5<pL5y + 5yD5<p^d p x (48) 
Again ignoring total divergences one finds 



m 1 d * 1+ y' 2 d 

M = — -— r — — — , 



r 2 dr Q 3 dr 

1 d 9 1 — d>' 2 d 
N = r , 

r 2 dr Q 3 dr ' 

1 d 2 0V d 



L = — — r 



r 2 dr Q 3 dr 

(49) 
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with L = I}. We can now rewrite 5 L as 

5 2 L = ^J d 3 x(50, Sy)H 



6y 



where 



and 



H=( M L 1=11^ 
\ iJ N J r 2 dr dr 



fr- 1 -^ " 1 + 0,2 ^ ^ deth-- 
~ Q { y'4>' i + y' 2 )' Q 4 ' 

The small fluctuation equation therefore becomes 

ii-^ = —z-r-r n—ip = uoip 
r 2 dr dr 

Again we first explore the existence of a zero mode ipo with 



with 



2a d i I a 
rh d-r^=[(3 



where a and (3 are constants. Setting 

and evaluating x/jq for the solutions of eq . (|46D we obtain with 

f°° 3 
ip = — I (J)' (x)dx 



the relation 

— ) I 1 . — , , , , , . 

-p / e V a 



e 

In the BPS limit with y' = 0' = £7 C , (5 = 1, the operator H of eq. 

H=-— r 2 ( - l ~ E c E - \ — 
r 2 dr \ E 2 1 — E 2 J dr 

Setting 

for an arbitrary function p(x) we have 

rr, I d 2 ,( -1-E 2 E 2 C 1 d 2 , / 

^ s "^ r M £ 2 i-£ 2 i "^ r M 
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Thus for arbitrary p(x), we have ip s Hip s = implying 5 2 L = or L constant in 
a specific direction about the BPS configuration. This behaviour may be inter- 
preted as indicative of a local symmetry, in this case of supersymmetry, and so of 
the cancellation of fermionic and bosonic contributions in the one loop approx- 
imation. Here, of course, we have no fermionic contributions and consequently 
those of the two bosonic fields have opposite signs. 



5 Fluctuations about the D— brane 



In the following we distinguish clearly between two different types of fluctuations. 
We consider the above BPS solution for the string attached to the 3-brane as 
background and consider first a scalar field propagating in a direction along the 
string and perpendicular to the brane and its anti-brane. The linearised equation 
of small fluctuations about this background is obtained from the second varia- 
tional derivative of the action which is the standard procedure and we therefore 

Our treatment here is somewhat different (see 
The resulting fluctuation equation has also been 



consider this first (cf. also [25 
below) from that in refs. 
given in ref.|2[]. It is necessary to return to the fully time-dependent version, i.e. 
1 



S 



x 



1 



(2tt)^ 

where in 3 + 1 dimensions 
only one scalar field y we have 
= d^y for i — 1, 2, 3 and fi = 



det(?7 Mi , + 



(59) 



In the electrostatic case with 



= (A , Ax, A 2 , A 3 , y, 0, 0, 0, 0, 0), F Qi 
1,2,3. Then 



Ej and 



det(?7^ + F f 



flU J 



-1 


Ex 


E2 


E 3 


d y 


-E x 


1 








dxy 


—E 2 





1 





d 2 y 


—E 3 








1 


d 3 y 


-d y 


-dxy 


-d 2 y ■ 


-d 3 y 


1 



(60) 



and so 



det (77 



[IV J 



:i-E 2 )(l + Vy 2 )-(B-Vy) 2 + (d oy y 



(61) 



We consider first the Lagrangian density (remembering that the relevant fields 
are A , Ai and y) 



C = l-Q, Q 



(1-E 2 )(l + Vy 2 ) + (E- Vy) 2 



y 2 



(62) 



The equations of the static Blon and the static catenoid discussed above follow 
again from the first variations 



dC 

dEi 
dC 

ddiy 
dC_ 

dd y 



1 

Q 



1 

Q 
1 

Q 



E l (l + Vy 2 )-d l y(V-Vy)_ 
^(l-E^ + ^E-Vy) 



(63) 
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In the BPS background given by 



diy = E t 



(Vyf = E • Vy 



El = yl Q 



(64) 



one finds 

8 2 C 
dEidE, 



8 2 C 



ddiyddjy 



8 2 C 



dEiddjy 



-EH. 



8 2 C 

8y 2 



(65) 

This enables us to write (ignoring again total divergences in shifting derivatives) 

8 2 C = (8A ,8A„8y)- 

( -diil + E^di 8,(1 + E 2 c )d 
8 (l + E 2 c )d, -d (l + E 2 )d 

+diEld Q 8,(1 - E 2 )di - d Q d J 



-diE 2 c di 
d E 2 di 



\ 



\ 



( SA 

SAi 



(66) 



In the linear approximation the Euler-Lagrange equations of the fluctuations 
5y = rj, 5Ei = 8o8Ai — di5A are therefore given by the following set of three 
equations 



-^ + d i (l-E%)d {n + d i & o (d 8A i -d i 8A< i ) = 0, 

j t (l + E 2 )(d 5A i -d i 5Ao)-j t E 2 c d i rj = 0, 
di(l + E 2 c )(d 8A i -d i 8Ao)- diE 2 c d iV = 



(67) 

(68) 
(69) 



The last of these three equations can be seen to be a constraint by appying d/dt 
and using the second equation. Substituting from the last 

diE^doSAi - dMo) = diEldiff - d^doSAi - $6 A ) 

into the first equation we obtain 



-ffjj + Ari- 8,(806 A ~ 8,5 A ) = 
The second of the three equations can be written in the form 

(1 + E 2 )(d 8A ~ 8i5A ) - E 2 c 8 iV = (1 + E 2 c )Q{r) 



(70) 



(71) 



where C(r) is an arbitrary function. Dividing eq. (|6Q]) by (1 + E 2 ) and taking the 
derivative 8,, we obtain 



8i(doSAi - 8i5Ao) = 8, 



El 



l + E 2 



diTj + diCi 



E 2 A A 2E C E' C Xi 

; Ar] + —r-^—diT] + did 



[l + El) 2 r 



(72) 



Replacing on the right hand side E 2 8ii] by the expression in eq.(|60"Dthis becomes 



8,(806 A, - 8,8 A ) 



El 



l + E? 



Ar] + 



2e:. 



E c (l + El) r 



[(8 5A, - 8,8 A c 



a] + did 
(73) 



15 



Choosing as gauge fixing condition the relation 



2E' 



[{doSAi - diSAo) - d] + did = 



E c (l + El) r 
one obtains the following fluctuation equation for r\ 







(74) 



All the relations from ( pQ) to (|74|) describe perturbations along the string and 
perpendicular to the brane. Eq. (|74| ) cannot be considered independently of the 
others as is apparent from the linkage of the fields in the above equations. Thus 
if one wants to determine the radiation of the string between the brane and the 
antibrane, one must connect the asymptotic behaviour of the field 77 with that of 
the vector field SA^. 

However, an equation like ([74]) is also obtained if one evaluates the determinant 
in the Born-Infeld Lagrangian at the BPS background and with an additional 
time-dependent scalar 77, representing the fluctuation field along a new spatial 
direction (cf. also ref.||). In this case this new scalar field in the D-brane 
background has no relevance to the string radiation, and we have 



det(?7^ + F f 



^v)\BPS,ri 



-1 


E x 


E 2 


E, 





d r] 


-E x 


1 








E x 


dif] 


—E2 





1 





E 2 


d 2 rj 


—E 3 








1 


E 3 


d 3 y 





-E x 


—E2 


-E 3 


1 





-do V 


-0177 


-d 2 rj 


-d 3 r) 





1 



(75) 



and so 



det(?7 All/ t ± liv)\BPS,rj 

Thus the Lagrangian density becomes 



-(i + E c 2 )(a 77) 2 -(^) 2 -i 



(76) 



(77) 



By expanding the square root and retaining only the lowest order terms, we again 
obtain a fluctuation equation like (|65|) , but this time for 77 with no relevance to 
radiation of the string. This is equivalent to studying the scattering of the scalar 
77 off a corresponding supergravity background. 



6 Absorption of scalar in background of D3 brane 

We now consider the equation of small fluctuations, i.e. eq.(|74]), in more detail. 
The fluctuation rj(t, x) represents a scalar field that impinges on the brane which 
reflects part of it and absorbs part of it depending on the energy u of the field. The 
absorption results from and takes place into the singularity of the real potential 
which corresponds to the black hole with zero event horizon in the analogous 
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case of the dilaton-axion system of e.g. ref. ||19|| . This absorption is a classical 
phenomenon. We therefore consider the equation 



(78) 



One can argue that the absorption is a consequence of the nonhermiticity of the 
potential. 

The radial part of this equation is with £ = r~ 1 ^>Yi m and angular momentum 



d 2 ^ 

dr 2 



+ 



1(1 + 1) 2 

• 5 + UJ 2 



= 



(79) 



This equation is a radial Schrodinger equation for an attractive singular potential 
oc r~ 4 but depends only on the single coupling parameter k = eu 2 for constant 
positive Schrodinger energy, i.e. for S'-waves the equation is with x = ur simply 



— I — ] o 

dx 2 x 4 / 



10) 



In the following we consider the general case, i.e. / ^ 0. The simplified case of 
the singular potential replaced by an effective delta-function potential has been 
considered in refs.0 and The solutions and properties of such equations have 
been studied in detail in the literature, in both the small- and large-/? domains 
and with inclusion of the centrifugal term — 1(1+1) / r 2 in eq. ( [75| ) for the calculation 

p7|| , ||28|| , [p9| . A recent investigation which 



of Regge trajectories I 
attempts to treat arbitrary power singular potentials is ref.[]30||. Eq. (|79|) describes 
waves above the singular potential well. With the substitutions 



vj/(r) = Ai/;(r), r = ^/ee z , h 2 = eu 2 , a = l + -. 



the equation becomes the modified Mathieu equation 

d 2 tfj 



dz 2 



+ 



2h cosh 2z — a z 



ip = 



51) 



12) 



which has been studied in detail in the literature |3]J (though some properties, 
such as large-h asymptotic expansions of Fourier coefficients, have even now not 
yet been published). Here we study the S-matrix in the domain of finite values 
of angular momentum / and h 2 ^ 0, i.e. in the domain of h 2 large. Relevant 



solutions and matching conditions for this case have been developed in [|32f and 
We follow the latter of these references here since this makes full use of 



the symmetries of the solutions. Moreover we can determine also the Floquet 



exponent v which ref. |32| leaves undetermined and only remarks that the notion 
that this is a known function of (our) a 2 and h 2 is "partly a convenient fiction" . 
For convenience we set in eq. fl82|) as in ref. f33| , |34fl 



-2h 2 + 2hq + 



A(q,h) 
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where q is a parameter to be determined as the solution of this equation and A/8 
is the remainder of the large-h asymptotic expansion (jS3| ), the various terms of 
which are determined concurrently with corresponding iteration contributions of 



the solutions ip of the equation and are known explicitly to many orders f34 
Then setting in eq. (|82|) 

ip(q, h; z) = A(q, h; z)exp[±2hi sinh z] (84) 

we obtain an equation for A which can be written 



dz 2 Ahi 



A d 2 A 
J ~ Iz 2 



i5) 



We let A q (z) be the solution of this equation when the right hand side is replaced 
by zero (i.e. in the limit h — > oo). Then one finds easily 

AJz) = -=L= ( 1 + i ^ z ) T<1 ' 4 ~» ^ e - z /2 e T*^ (g6) 
Vcoshz V 1 - «sinhz/ 

Correspondingly the various solutions tp are 

ip(q,h;z) = A g (z)e^[±2/Msinhz] exp ^ ±lhe ^ e ^/^ 

V cosh z 

= A,(z)esp[±2tasinhz] exp (T ^ e ' Z |) e ^/ 4 (87) 

vcosh^ 

We make the important observation that given one solution ip(q,h;z) we can 
obtain the linearly independent one either as ip(—q, —h; z) or as ip(q, h; —z), the 
expression (|83|) remaining unchanged. With the solutions as they stand, of course 
ip(q,h;z) = ip(—q,—h;—z). Below we require solutions He^(z),i = 1,2,3,4, 
with some specific asymptotic behaviour. We define these in terms of the function 

Ke(q, h; z) := ex M CT g/ 4 ] A q {z)exp[2hi sinh z] = k{q, h)ip{q, h; z) (88) 
V—2ih 

Since this function differs from a solution ip by a factor k(q, h), it is still a solution 
but not with the symmetry property ip(q, h; z) = ip(—q, —h; —z). Instead, after 
performing this cycle of replacements the function picks up a factor, i.e. 



in leading order. One can easily show that the quantity $o of ref. [g2| is related 
to q by $o = iqn/2 + 0(l/h). In Fig. 1 we show the behaviour of q as a function 
of h. In order to be able to obtain the S— matrix, we have to match a solution 
valid at z = — oo to a combination of solutions valid at z = oo. This is achieved 
with the help of Floquet solutions Me± u (z, h 2 ). As such, these satisfy the same 
circuit relation as a solution M±l(z, h 2 ) of eq.(|82|) expanded in a series of Bessel 
functions, i.e. we have the proportionality 

Me u (z, h 2 ) = a v Ml 1] {z, h 2 ), a u (h 2 ) = Me u (0, /i 2 )/MW(0, h 2 ) (90) 
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The functions Me± v (z, h 2 ) are expansions of the modified (hence 'M' instead 
of 'm') Mathieu equation in terms of exponentials (hence 'e') which are uni- 
formly convergent in any finite domain of z. For large values of the argument 
2h cosh z of the Bessel functions of the modified Mathieu function Mfp (z, h 2 ) can 
be reexpressed in terms of Hankel functions. With the dominant terms of these 
we can obtain the large 2/icoshz asymptotic behaviour of the Floquet function 
Me± u (z, h 2 ), i.e. for \z\ — > oo 

. , r , , cos(2/i cosh 2; =F vti/2 — 7r/4) , n . 

Me±„{z, h 2 ) ~ exp[±z7r 7 /2] ^ ' '-J. (91) 

V In cosh z 

where (with Me u (-z,h 2 ) = Me. u (z,h 2 )) 

expfrwy] = = M«(0, h 2 )/M^(0, h 2 ) (92) 

We now define the following set of solutions of eq. fl82|) by setting 

He {2) (z,q,h) = Ke(q,h,z),He w (z,q,h) = He {2 \z, -q, -h), 
He (3 \z,q,h) = He m {-z,q,h),He^{z,q,h) = He {2 \-z,q,h) (93) 

The solutions so defined have the following asymptotic behaviour (where e(z) = 
(2hcoshz)- 1 / 2 ): 

He^(z,q,h) = e(z) ■ exp[-ihe z - i^}, 'Rz » 0, 
r _*oo exp[— iur — m /A\ 



7T. 



He i - 2 \z,q,h) = e(z) ■ exp[ihe z + i—], Uz » 0, 
r _,.oo exp[iur + in / 4] 



He^ 3 \z,q,h) = e(z) ■ exp[-ihe^z\ - i-}, << 0, 

He^\z,q,h) = e{z) ■ exp[ihe^z\ + i-], « 0, 
r ^o r x l 2 exp\ieu) jr + z~] 



(e^)V2 



(94) 



For the following reasons we choose the latter, i.e. the solution He^ 4 '(z, q, h), as 
our solution at r = 0. The time-dependent wave function with this asymptotic 
behaviour is proportional to 

g— iut+ieui/r+iir /4 

Fixing the wave front by setting (p = —ut + euj/r + n/A = const, and considering 
the propagation of this wave front, we have 

eu 

T = 

(f + Lot — 7r/4 
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so that when t — > oo : r — > 0. This means that the origin of coordinates acts as 
a sink. 

With eq. ( |9T| ) we therefore equate in the domain ^kz » 0: 



MeJz,h 2 



Me_ u (z, h 2 ) 



-exp[iiT'-y/2} ^expliv—jHe^^z, q, h) — exp[—iv^]He^ 2 \z, q, h) 



-27T7/2] jexp[— iv— ]He^\z, q, h) 



-exp[ 
exp[iu—]He^ 2 \z, q, h) 



(95) 



where the second relation was obtained by changing the sign on v in the first. 
Changing the sign of z we obtain in the domain 9ftz << 0: 

Me u (-z,h 2 ) = Me- U (z,h 2 ) 

— -exp[i7T~f/2]<exp[iv—]He ( - 3 \z,q,h) — exp[—iu—]He^\z,q,h) 



Me v {z,h 2 ) 



-exp[—i7cy/2]lexp[—ib'—]He( 3 \z, q, h) 



7T 

exp[iu— }He^ (z, q, h) 



(96) 



These relations are now valid over the entire range of z. Substituting eqs.(|96|)into 
eqs.(p5|) and eliminating He^ we obtain 

— Sin 7TZA = sm7r ( 1 + is).He w (z,q,h) — smn^.He^ 2 \z,q,h) (97) 
In a similar way one obtains the relations 

— sin iris. He( 2 \z, q, h) = sin7r in + u).He^\z,q } h) — sin Ti'j. 
smirv.He^(z,q,h) = —sin tcj. He (3 \z,q,h) + sin7r(7 — v).He^ A \z, q, h) 

(98) 



^From eqs. (|89|) and (|9~3D we see that He^(z, q, h) is proportional to He^ 3 '(z, q, h). 
^From ( p9|) and ([58") we see that the proportionality factor is given by 



exp[i^{q + 1)] = _ sm7 [( r ) / + zy ) 
2 sm ixv 



(99) 



^From eq. fl97]) we can now deduce the S-matrix Si = e 2lSl , where 5i is the phase 
shift. The latter is defined by the following large r behaviour of the solution 
chosen at r = 0, which in our case is the solution He^\ Thus here the S-matrix 
is defined by (using fl9~7|)) 



smniy- 



^,l/2gieo;/r+i7r/4 

(ecu) 1 / 2 



sin7r(7 + u)e J7r / 4 


"sin7T7(— 1)' 


\Jurr 


.sin 7r(7 + v) 



-e in / 2 e iujr 



:-ire- 



e e 



iS p—iln/2 



2i^¥ 
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(100) 



^From this we deduce that 



Si 



sm7r T c wr(t+l/a) 



sin 71(7 + v) 



smnu 



(101) 



We can see the relation of this high-energy (i.e. large \h\) expression of the S- 
matrix to the low-energy expression of ref.|U[ by recalling that R of the latter 
is here exp{m r y). With this identification we can write S\ 



Si 



R ~7i 



^iir(l+l/2) 



(Re 1 ™ - 



R 



which agrees with the S'-matrix of ref.fJHJ, i.e. we thus obtained the same exact 
expression of the S-matrix here with our large-/i considerations. In fact, compar- 
ison with the considerations given there allows one to write down the reflection 
and transmission amplitudes A r and A t as A r = 2i sin 717 and A t = sin ttu re- 
spectively. We thus have one and the same expression for the S"-matrix for the 
two asymptotic regions, i.e. in the low energy and high energy domains. One 
should therefore be able to proceed directly to the large-/t case from the exact 
S'-matrix derived in the small-/i domain. This is an interesting calculation which 
we do not attempt to go into here. We only indicate in Appendix A the first 
necessary step in that direction, i.e. the derivation of large-A asymptotic ex- 
pansions for the Fourier coefficients of Mathieu functions. In this connection we 
make the following two observations. 1) Eq.fl80"D is invariant under interchanges 
x <-> k/x,^ <-» ^>x which means that the inner or string region is equivalent 
or dual to the outer or brane region. 2) Due to the SL(2,R) invariance of the 
.D3-brane its action is mapped into that of an equivalent Z)3-brane by S— duality 



transformations [35| or, in other words, weak-strong duality takes the D3-brane 
into itself ||. It would be interesting to find some connection between these 
properties, or equivalently the symmetry which the SL(2, R) invariance of the 
D3-brane action imposes on the S— matrix. 

The quantity 7 is now to be determined from eq. (p9]). One finds 



sin 7T7 = sin nu 



ie 1 " 29 cos ttu ± 



1 + e tnq sin iris 



(102) 



It remains to determine the Floquet exponent v in terms of q and h. In Appendix 
B we derive the appropriate large-/i behaviour of v for the case of the periodic 
Mathieu equation. Replacing there the eigenvalue A by a = (1+ 1) 2 and observing 
that h 2 remains h 2 , the appropriate relation for our considerations is 



cos iris + 1 



7lC 



Hi 



{8h)i/ 2 
(8/1)9/2 



1 + 



3(g 2 +l) 
64h 



n 



1+ 



4J u 4 

3(" 2 

Uh 2 



+ 0( 



h 2 ' 



g+i))r(4 i )cos(^^ 



2tt29/ 2 



(103) 



Since the right hand side grows exponentially with increasing h the Floquet ex- 
ponent v must have a large imaginary part. Since the right hand side is real, the 
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real part of v must be an integer. Using Stirling's formula we can approximate 
the equation for q ~ h (i.e. irrespective of what the value of I is) as 



cos7rz/+ 1 = J~ cos(^)(e732) /i/2 
^From eq.( p.01| ) and eq. (|102|) we obtain 



^cos(^) (104) 



Si = ie dn ( cos ttu — v cos 2 nu — 1 — e iq7T 



(105) 

^From this we obtain the absorptivity A(l, h) of the 1-th partial wave, i.e. 

A(l,h) : = 1 - |^| 2 (106) 



with near asymptotic behaviour 



A(l,h) ~ 1 



2tt(16/i) 9 



,8/i 



(107) 



In Figs. 2, 3 and 4 we plot A(l, h) as a function of h. One can clearly see the 
expected asymptotic approach to unity and in Fig. 2 some sign of rapidly damped 
oscillations. This behaviour agrees with that obtained on general grounds in 



ref. ||15|| . We also observe that in the high energy limit logarithmic contributions 
as in the low energy expansions, discovered originally in 17, 18], and typical 



of the low energy expansions of |L3j and jnj, do not arise. Of course, these plots 
do not extend down to h = 0, since our asymptotic solutions become meaningless 
in that domain. The continuation to h = can be obtained, however, from small- 



h expansions such as those derived in refs. fl~3H and |E{|. Thus the absorptivity 
A(l, h) is known over the entire range of h. We observe that Si = for q = 
1, 3, 5, • • •, with [(/ + 1/2) 2 + 2h 2 }/2h ~ 1, 3, 5, • • •. Only in the plot for I = 2 is h 
sufficiently large to hint at these zeros. 



7 Concluding remarks 

Branes, whether fundamental or solitonic, play an important role in all aspects 
of string theory. In particular D-branes have been looked at as string-theory 
analogues of solitons of simple field theories, and some of their important proper- 
ties such as charges are well understood. Our first objective in the above was to 
investigate properties of solitonic objects of Born-Infeld theory in ways familiar 
from field theory, in particular their classical stability. It was shown that the Blon 
and the catenoid as distinct, i.e. free objects, are stable configurations whereas 
the brane-antibrane system is unstable; we also recognised the zero modes as- 
sociated with these and their significance. We then considered the D3-brane of 
Born-Infeld theory and recognised this as a BPS state that preserves half of the 
number of supersymmetries as discussed in detail already in |fj] . The equation of 
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small fluctuations about this D3-brane was derived and shown to be convertible 
into a modified Mathieu equation. The low energy solutions of this equation, the 
S'-matrix for scattering of a massless scalar off the brane and the corresponding 
absorption and reflection amplitudes are similar to those for the dilaton-axion 
system investigated first in refs. fl6|, [17], HU, where the important logarithmic con- 
tributions were discovered, and then investigated in extensive detail in O] and 



19|j . Here we performed the high energy calculations which complement in par- 
ticular those of |19[, thus completing the investigation of the modified Mathieu 
equation for the purpose of obtaining absorption cross sections for all such cases. 
In particular the behaviour of the important Floquet exponent involved in these 
calculations (in general a complex quantity) is now fully understood, the Floquet 
exponent being vital in the evaluation of the S-matrix which we derive and the 
calculation of the corresponding absorption amplitudes and cross sections. Ac- 
cording to our findings the high energy limit of the absorption cross section does 
not involve logarithmic contributions, quite contrary to the low energy limit. 

The high energy case considered here is not only of interest in the immediate 
context of the Born-Infeld model considered here, but together with the low- 
energy case also of considerable interest in connection with the concept of duality 
which links weak coupling with strong coupling. The D3-brane with Schrodinger 
potential coupling eu 2 , which links the gauge field charge e with energy^ of the 
incoming scalar field is presumably the ideal example for the investigation of 
this property. Investigations elucidating this aspect are of considerable inter- 
est. We also envisage interest in the study of non-BPS configurations, including 
sphalerons and bounces, as a matter of principle, i.e. even if the effect of these 
is not of dominant importance. Finally we remark that it should be possible to 
proceed directly from the S'-matrix derived in ref.|]TU| to the high-energy case 
here by using appropriate asymptotic expansions for the cylindrical functions 
and expansion coefficients involved (for the latter such expansions do not seem 
to have been given in the published literature so far, but we comment on these 
in Appendix A). 
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Appendix A 



In ref. |19| on the absorptivity of the L>3-brane of the dilaton-axion system 
it was shown that the S-matrix for scattering of a massless scalar field off the 
brane is given by 



S 



(R - i) e -*rf 



Re 1 



(A.l) 



R 



where 



R 



Mi 1] (0,h) 



Mjp(z, h) being the modified Mathieu function expanded in terms of Bessel func- 
tions, i.e. 

oo 

Me v {0,h)M^{z,h)= ]T c u 2r (h 2 )J u+2r (2hcoshz) 

r=— oo 

(an expansion with better convergence to use in practice is one in terms of prod- 
ucts of Bessel functions as shown in ref. |19|]) where Me y (z,h) is the Fourier or 
Floquet solution of the Mathieu equation. In the published literature the co- 
efficients c 2r (h 2 ) have only been considered as power series in rising powers of 
h 2 , and consequently were used in ref. 
main. 

high energy case directly from this expression by developing large-/; 2 asymptotic 
expansions of the Mathieu function Fourier coefficients c 2r (h 2 ) (for the Bessel 
functions the corresponding expansions are known). We know of no publication 
where such expansions have been given, but one of us (M.-K.) remembers from 



19 



in the small h 2 or low energy do- 
It would be very interesting to make the transition to the large-/i 2 or 



private communication with the author of ref. |37J that these Stokes-type asymp- 
totic expansions can indeed be obtained. One writes the recurrence relation of 
the coefficients (cf. [EH]], p. 106) 



[A - (v + 2p) 2 } 

C2p+2 + £2p-2 — ?2 c 2p 

(the Mathieu equation being y" + (A — 2h 2 cos2x)y = 0). For \h 2 \ 
implies 

•(2p+2)/2 



(A.2) 
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Setting 
we have 
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^From this we deduce that the next approximation to C2 P +2 is obtained from 
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The sums on the right hand side can be evaluated. E.g. 

o 1 



Ep 2 

p=o 

so that one obtains 



l 2 + 2 2 + 3 2 
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r(r + l)(2r + 1) 



Pr 
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r(r + l)(2r + 1) ^ r ( r 



6 
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(A.5) 



Proceeding in this way one can indeed obtain the desired asymptotic expansion 
of the coefficients c<i p - (In fact the asymptotic expansion of the Bessel function 
- similar to that of a linear combination of Hankel functions - can be obtained 
from its recurrence relation in a very similar way). 



Appendix B 

For the determination of the large-h behaviour of the Floquet exponent v we 
make use of results of ref.[|34|]. A fundamental pair y I: y u of respectively even and 
odd solutions of the original periodic Mathieu equation with eigenvalue A defined 
by 

Vi{-z) = yi(z), y n (z) = -y n (-z) 



can be chosen to satisfy the following boundary conditions (cf. e.g. plfl , pp.99, 100) 

yi (0) = l, y/7(0) = 0, y/(0) = 0, y//(0) = l 

^From its original defining property the Floquet exponent v can then be shown 
to be given by (cf. pi, p. 101) 



cos tcu = yi(ir, A, h 2 ) 



so that (cf.[ 



p. 100) 

COS TTV + 1 



2y I (7i/2;X,h 2 )y I /(7c/2;\,h 2 



(b.i; 



(B.2) 



The solutions yi(z),yn(z) can be identified with the large-h solutions ce, se of 
ref . p4l (there eqs.(64))in terms of functions A(z),A(z) as in eq.fl84D above with 
normalization constants N , Nq, i.e. in leading order 

ce(o) = 2A A(0), se'(0) = AhN 'A(0) 

from which we deduce in leading order for large \h\ that 

A = 2" 3 / 2 , N ' = 2-^ 2 /h 

Eqs.(65) of ref.|34]] give the large-h expansions of ?//(7r/2; A, h 2 ) and yj/(n/2; A, h 2 ). 
Inserting these multiplied by the appropriate normalization constants into eq.( |B.2|) 
and retaining the dominant terms for large \h\ we obtain 



cos tip + 1 



{h 



]h)i/ 2 



(B.3) 
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in agreement with a result cited in ref. |TT|(p.210) from with logarithmic cor- 
rections. We, however, see no such logarithmic terms in the simpler formulation 
of ref. [[54]]. The relation ( |B.3| ) we rediscovered here has practically been unknown, 
largely in view of the difficulty to extract it from the complicated considerations 
of ref.f36|. Our derivation above is simple and closes a difficult gap which the 
author of ref. [ 32 1 commented upon with the words: "It is not likely at this stage 
that an analytic relation will ever be found connecting (our) v and 7 to (our) a 2 
and h 2n . Our search of later literature did not uncover other derivations. The 
main source summarizing more recent developments in the field of the Mathieu 
equation is ref. [[J^ . 
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Figure Captions 



Figure 1 

The function q(h) plotted versus h, which, of course, is valid only away from 
h = 0. The plot should be compared with graphs in ref.[[J2|] where a similar but 
less convenient quantity is used. 

Figure 2 

The absorptivity A(l, h) for I = 0. 

Figure 3 

The absorptivity A(l, h) for I = 1. 

Figure 4 

The absorptivity A(l, h) for I = 2. 
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